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Abstract 

A rigorous and most general covariant kinetic formalism is developed to study the nonlinear waves 
interaction in relativistic Vlasov plasmas. The typical nonlinear plasma reaction is a nonlinear 
current measured by the nonlinear plasma conductivity, and these quantities are derived here on 
the basis of relativistic Vlasov-Maxwell equations. Knowing the nonlinear plasma conductivity 
allows us to determine all plasma modes nonlinearly excited in plasma. The general covariant form 
of nonlinear conductivity is provided first for any value of plasma temperature and for the whole 
complex frequency plane by a correct analytical continuation. The further analysis is restricted to 
a correct relativistic particle distribution which is vanishing for particle speeds greater than speed 
of light. In the limit of nonrelativistic plasma temperatures the covariant nonlinear conductivity 
is significantly different from the standard noncovariant nonrelativistic results which are reached 
only in the formal limit of an infinitely large speed of light c — > oo. 
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I. INTRODUCTION 



Under various sources of free energy plasma becomes a high nonlinear dispersive medium. 
Numerous nonlinear wave processes occur in both laboratory and space plasmas [Jlj]: e.g. in 
laser plasma interaction the high intensity pump waves stimulate the parametric instabilities 
and the nonlinear collisionless damping of plasma waves; or the shock structures arising in 
many astrophysical plasmas fed the well-known two-stream (electrostatic or Weibel) insta- 
bilities [Q]. 

Even for relativistic temperatures the collisionless (Vlasov) plasma condition applies to 
many astrophysical scenarios where dissipation is dominated by wave-particle interactions 
rather than binary collisions. In these cases a fundamental kinetic description is required, but 
the existing standard nonrelativistic results could be improved by the covariant relativistic 
approaches which reproduce the general and relativistically correct dispersion relations in 
any frame that is not necessarily inertial. 

In the first earlier treatment [|3|] (named as Paper I in the next) the nonlinear three-waves 
interaction have been investigated in the case of relativistic Vlasov plasma with isotropic 
distribution function of particles. Using the relativistic Vlasov-Maxwell equations the new 
relativistic forms of nonlinear current and nonlinear plasma conductivity have been derived. 
The nonlinear current is the typical nonlinear response of a plasma medium where a pump 
and an idler wave interact exciting a third signal wave. 

The nonlinear current seems to play an important role also in the saturation stage of 
plasma filamentation arising in laser plasma interaction [0|] or in the jet-plasma astro- 
physical structures [Js]]: shocks in GRBs sources, pulsar outflows, or solar winds. The 
filamentation instability is the final stage of nonlinear waves interaction when the plasma 
distribution function is anisotropic and which persists as long as the external perturbation 
is sufficiently strong. This case should be analyzed using an anisotropic bi-Maxwellian dis- 
tribution function to calculate the nonlinear conductivity, and it will be discused in the next 
papers of this series. 

In the present paper we continue to develop a fully relativistic and covariant kinetic 
theory of nonlinear plasma waves interaction providing an analytical model to estimate the 
plasma nonlinear response. Starting from the Maxwell equations we show in Section II that 
two waves coupling leads to a nonlinear component of plasma current which generates the 
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third wave signal. To find the corresponding nonlinear fluctuation of plasma density we use 
as in Paper I the relativistic nonlinear Vlasov equation to obtain the general form of the 
nonlinear current arising from the interaction of two waves with any polarization and any 
direction of propagation. 

Considering only longitudinal waves the general covariant expression of nonlinear plasma 
conductivity is derived in Section III for any subluminal or superluminal phase velocity 
and for the whole complex frequency plane. The nonlinear conductivity is then calculated 
in Section IV using an appropriate relativistic distribution function which is vanishing for 
particle speeds greater than speed of light. Knowing the nonlinear conductivity we could 
find the electric fields of interacting waves as solutions of nonlinear (coupled) equations 
system (42)-(44) from Paper I. 

For nonrelativistic plasma temperatures we derive in Section V apparently for the first 
time, the covariant expression of nonlinear conductivity in terms of well documented plasma 
dispersion function We show that the covariant form is markedly different from the 

standard classic expression of nonlinear conductivity provided by the noncovariant nonrel- 
ativistic theory and which can be obtained only in the unphysical limit of infinite speed of 
light c — > oo (because the classical theory assumes all waves as being subluminal). 

II. BASIC FRAMEWORK 

The basic model followed in this paper is an infinitely extended collisionless plasma of 
high temperature electrons which requires that the electrons be treated relativistically. For 
plasmas with low collisionality [y c juj v ~ g << 1), the cooperative motion is due to the 
electromagnetic coupling of the particles, not to collisions. Therefore, the relativistic Vlasov- 
Maxwell equations are used, and for plasma fluctuations with sufficiently large amplitudes, 
the nonlinear forms of them are required (see notations and equations (13)-(17) from Paper 



which follows directly from the Maxwell equations (see equations (14)- (17) from [3]). Waves 



I). 



The equation governing wave propagation in plasmas, is: 




(1) 
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interaction gives rise to the nonlinear terms in J: 

J(r,t) = J L M)+J NL (r,t), (2) 



with 



J L (r,t) = 3iQk,"i) = (3) 



and 



J NL (r, t) = J2 J? L ( r . t) = E J « L ( k - ^) = E JT^ M - ( 4 ) 

p>2 « « 

We assume that 

E(r, t) = E E *( k - <«*) = E E°e-^" k - r \ (5) 

i i 

D(r, t) = eE(r, t) = E e (^) E *( k *. ( 6 ) 

i 

and E° is taken as essentially independent of position in space. 

For the lowest-order nonlinear phenomena (p = 2), as typical three- waves configuration, 
two waves (indices 1 and 2) give rise to the third one (index 3) through a first generated 
nonlinear current: 

jf(r,t)=jf 2 L (r,t) + J^(r,t). (7) 

With ©-© and 



e(ui) = 1 + a(ui), (8) 



(JTJ) becomes: 



1 d 2 

Vx(Vx) +e -- 



2 1 ^ flTNL 



clearly indicating (with Q and 0) that wave electric fields E*(ki, UJi) are coupled through 
nonlinear current J NL . 

A. Nonlinear current 

The nonlinear density fluctuation / NL is solution of the nonlinear Vlasov equation (see 
in Ref. [ 311) and it will determine the nonlinear current: 
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T NL 
J 3 



a NL EE = q d 3 pvf NL 

J — oo 

p 



— iq 



oo 

2 / j3 



d A p- 



Ei + PXBA% 
7mc / ap 



(10) 



7m^ 3 - k 3 • p 

where (1 <-> 2) denotes the interchanging of indices 1 and 2 required by (J7J), and 7 = 
y/l + p 2 /m 2 c 2 is the Lorentz factor. Otherwise, the linear fluctuations of the distribution 
function are solutions of the linear Vlasov equation as 



fi = -i 



Q 



dfo = q dfp p Ej 
<9p <9p p 



(11) 



where an isotropic distribution function /o(p) = fo(p), has been assumed at equilibrium, 
and the relativistic resonance factors are given by 



Ri = 1 - (p • ki/TmUi) (i = 1,2,3). 



From (JTTJ) 



(12) 



<9/i • 9 ^/o fp-Ei 
— z- 



<9p oJiRiP dp y p 



p d 2 f /dp 2 | 1-^-7^ p + 



dfo/ dp 



I - Ri pki 



+ E, 



p ■ k 

and a new expression for density current (fTUj) is derived 



JNL = C r NL E 1 E 2 



2vrg 3 /- /"» pp ± a/ 
ap|i / ap_L — 



m 



oo Jo 



ft 



7P <9p uiu 2 uj 3 RiR2R3 



E^E 



2- 



with the second order tensor 



(13) 



ft 



V 



d 2 f /dp 2 
dfo/dp 



- J (WiiZj + W 2 #2 
;k!-k 2 )c 2 



tO\U)2 



I- 1- 



Ri 
pki 
RiJ P ■ ki 



1 



PP 

2 



1 

R~2 



k 2 p 
P • k 2 



(14) 
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Comparing with Eq. (25) from Paper I describing only the interaction of parallel propagating 
waves, kj || k, = 1, 2, 3), the new general expression (fHIj) (combined with (JHJ)), provides 
the nonlinear current arising from the interaction of two waves with any polarization and 
any directions of propagation. 

But for the sake of simplicity we consider in the next only the case of electrostatic 
(longitudinal) turbulence (Ej || kj, i = 1,2) with interaction of parallel propagating waves 
(along the same direction) ki || k2. 

III. CO VARIANT CASE 

In the general covariant case it is convenient to transform to the new variables of integra- 



tion [j7j], y = p\\/(mc) and E = y 1 + ^ 2 J; , which yields for the corresponding components 
of nonlinear conductivity tensor in (j!3j) 



a 



^._ Wm ^^§/_^ d ,yi^_^_. (16) 

where Qll is the component of second order tensor f2 corresponding to the interaction of 
two longitudinal waves (Ej || k i; i = 1,2): 

n LL = £ {uiRl+aA) + + (17) 

and 



E \dE 2 ' dEJ v * \ E 2 J \ R 2 Ri 



fl ' = 1 -S = 1 -i (i=w) ' (18) 

In (118)1 the inverse index of refraction z = u/kc is introduced. 



A. Reduction of aF L 
Introducing (|T7|l in (fTK|) we obtain 

<rf L = -27rg 3 mc 2 (/ 1 + J 2 ) (19) 
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with 



dEd% r^ 1 3 i 

1 " !, E 2 dE* J_ vw - 1 Vy u^R.R, 



+ (1 <->2) 



(20) 



dE df 



h " L E dE dyV y 1 E 2 J u^RlR, 



ir 



+ (!<-> 2) 



(21) 



where (1^2) denotes the interchanging of indices 1 and 2. In terms of the functions 
u = u(E) and v = v(E) defined by 



du 
dE 



-E 



,d 2 f dv 



dE 2 ' dE 'dE 



and using the substitution t = y/E the integrals (j2Uj) and (}2*T]) from above read 



(22) 



h 



^1^3 Jl 



dE 



du 



Vi-e- 



dt 



dE J_^E^ Oi - t) (z 3 - t) 



+ (1 - 2) 



2zi 



^1^3(2:3 -Zi) A 



dE 



Uy/W 



E 2 



1)E 2 + 1 



l)E 2 + 1 



+ (1 <-> 2) (23) 



z\z 3 



^1^3 Jl 



2zfz 3 



^1^3(^3 - z x ) 2 J 1 



dE 



dEj_ VT ^ F ? ( Zl - t ) 2 (z 3 -t) 
vVE 2 - 1 f 2zi(z 3 - ^)£ 2 



£ 4 [ [{zl - 1)E 2 + l] 2 {z\ - l)E 2 + 1 



From 



{zl ~ l)E 2 + 1 



+ (!<-> 2). 



(24) 



(25) 



which is used in (123)) to find 



/1 



2*i 



^1^3(2:3 - zi) e 2 Ve 2 - 1 



>f - 1)£ 2 + 1 (z 2 - 1)£ 2 + 1 



+ 



2 / dE 



jVE 2 - 1 



[(z 2 - l)E 2 + l] 2 [{z 2 - l)E 2 + l] 2 



+ (1^2) (26) 



and collecting ()26|) and (ffifj) the component ()19j) of nonlinear conductivity finally reads 



.jv L 4:7iq 3 mc 2 Z\ 



X 



2:3 - z x 



1 + 



^3 
2 3 - Zi 



dE 



E A 



4W - 1 



1 - zf)E 2 - 1 



dE 



>{E)VE 2 - 1 
E 2 



*3 



+ (!<-► 2) 



(27) 



.[(l-z 2 )£ 2 -l] 2 [(l-2 2 )£ 2 -l] 2 _ 
with v(E) defined by the isotropic distribution function fo{E) in (j22j). 

Relation (|2*Tj) is the general covariant form of longitudinal conductivity plasma response 
to the interaction of two (longitudinal) plasma waves. We can use it to find the expressions 
of plasma waves solving the nonlinear equation system (42)-(44) from Paper I. But relation 
()27|) is not valid for the whole complex frequency plane because it has been obtained with 
assumption of a positive imaginary part, T > 0, of frequency uj = uj t + iT (see also in Paper 
I). The index "+" in (j2"Tj) indicates this condition. 



B. Analytic continuation for subluminal waves 

As we noted, the relation (}2Tj) holds for positive values of the imaginary part of the 
frequency V = Quj > 0, corresponding to Qz > 0. In order to derive the corresponding 
dispersion relations for negative values of the imaginary part of frequency, T < (i.e. 
Qz < 0), we have to analytically continue the integrals in ()27|) into the negative imaginary 
plane of the complex variable z. 

For superluminal waves generally defined in unmagnetized plasma by \z r \ = \v p h ase /c\ > 1, 
all the integrals in (jTTj) admit no pole inside the integration interval, and we simply have 

~,NL „+,NL {r)Q \ 

which also holds for plasma waves with subluminal phase velocities \z r \ < 1 but with sus- 
ceptible growing amplitudes T > 0. 

For subluminal damped waves with \z r \ < 1 and T < each of the integrals in (f2*7j) 
admits the pole E ci = 1 / a/ 1 — zf > 1 inside the integration interval [1, 00], and the 
analytical continuation of (|2*Tj) into negative imaginary plane (negative index) reads as 
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follows 



— ,NL +,NL . 
a \\ = °V + 



4:TT 2 q 3 mc 2 Z\ 



U1U3 z 3 - Zi 



v(E c3 )zU 1 - Z; 



:i - zt) 1 + 



Z\Z 3 
Z3 - Zl 



v(E c i)zlJ 



l-zi 



1+ 



Z\Z 3 

z 3 - z x 



-l + *3 



+ (!<-> 2). 



(29) 



C. Interlude 



Combining now (|27)l and (|2*9"|) we derive the general form of the longitudinal conductivity 



NL +,NL . 4:7r 2 q 3 mc 2 z x 

On — (7,1 + 2 X 

11 11 a;ia;3 23 - Zi 



#(1 - i^Di^-^M^KV 1 - z 



1 - 1 1 



-#(1 - |9ffci|)ff(-92i) v(E cl )4Jl - z 2 



;i - *?) 1 + 



- «i 

^1^3 



^3 - Zl 



-1 + Z3 



+ (1 <-> 2). 



(30) 



holding for the whole complex frequency plane, i.e. for interaction of any growing or damped 
waves. The frequency and the wave-number of signal wave (noted here with index 3) are 
unknown in the conductivity expression, but they are usually provided by the conservation 
laws 



uj 3 = u\ ± lu 2 , k 3 = ki ± k 2 . 



(31) 



IV. RELATIVISTIC MAXWELLIAN PLASMA 



The isotropic equilibrium distribution function f is supposed to be Maxwellian. But 
a rigorous relativistic analysis requires for an appropriate relativistic distribution function 
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which is vanishing for particle speeds greater than speed of light. We therefore consider the 
Maxwell-Boltzmann-Jiittner ditribution function 



Up) = Ce 



-y,E 



(32) 



with a normalization constant 

C 



mc 



(33) 



so that f d 3 p /o = no- K u {fi) denotes the modified Bessel function. We use in the 
second equation of to find 



v = -(E + -)Ce~ tlE , 
A 4 



which can be introduced in (|2 



(34) 



a 



,NL 



Anq mc Z\ 



^1^3 z 3 - Z\ 



1 + ^^)J 1 -2J 2 

z 3 ~ Z\ 



+ (1^2) 



to calculate the integrals: 



(35) 



Jt 



~» E VE 2 -l / 1 



£ 4 



;i - zl)E 2 - 1 (1 - zf)E 2 - 1 



(36) 



4 = 

- Cj r«£^E!( i+ I) 

V. NONRELATIVISTIC THERMAL PLASMAS (yu >> 1) 

We now consider the limit of nonrelativistic plasma temperatures /i » 1, and to lowest 
order in /i" 1 << 1 

KM*<[^e-M C=(£)Wj*L s e>>. (38) 
V 2(/x) 27r (mc) 6 
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*3 



.[(l-sl^-lP [(1 1]* 



(37) 



and with the substitution E = a/1 + s 2 the first integral (|36j) becomes: 



2tt' 



"OO 

•'' ' ds 



mc 



g 2 e -nVi+s 2 



1 + s 



2\2 



(i - 4)# - 4 (i - #y - z\ 



(39) 



Because of the exponential function the main contribution to the integral for large values of 
fi a >> 1 comes from small values of s « 1, so that we may approximate a/1 + s 2 — 1 ~ s 2 /2 
yielding 



JJ_\ 3 / 2 n 
2tt, 



mc 



1-z 



„, .2 /-oo p -^ 2 /2 



\ - 7 2 , 

1 z 3 JO 



l-zi 



(3-1)} 



Uq/jL 



2 3 / 2 7r(mc) 3 



f 3 Z+(f 3 ) hZ+{h 



i-4 



+ 



Z-x Z-i 



(40) 



1-2? (1-^ 2 )(1-% 2 )J 
where we have used (in the positive imaginary frequency plane) the plasma dispersion func 

tion of Fried and Conte [p] 



7T 



-1/2 



dx -. 

oo X J 



(41) 



of the argument 



fi 



2^ 



And the second integral (|37j) is calculated in the same manner as 
n /i 



J- 



2 3 / 2 7r(mc) 3 [ 4(1 - z 2 ) 



hZ+{h/2) (1 - ^ 3 )A^ + (/i/2) + z 2 -z 2 + z 3 (l - zl 



4(1-^; 



2(1-^) (l-*g 



(42) 



(43) 



For nonrelativistic plasma temperatures > 1) we are entitled to assume |/| ^> 1 for 
which (see the asymptotic approximation of plasma dispersion function [J(|)] Z + (f/2) ~ 
2Z + (f ) and after summation of the integrals ()4U}) and (J4*3"j) in ()35|) we obtain 



+,nl n q 3 /j f z x z % 

On — 



V2m 2 cuiuj 3 z 3 - zi V 23 - Zi 



1 + J + + (1 - 2) (44) 
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with J + given by 



J+ = / 3 Z+(/ 3 ) _ (l + z 3 )f 1 Z+(f 1 ) + z\ - z[ - z 3 (l - z 



1-4 



i-z\ 1 (1-zfKi-zi) ■ (45) 

For negative imaginary frequencies we can find cr^ ,NL substituting (|H4|) in (j29jl . or it is 
easier to use in (|45j) the well-known analytic continuation of the plasma dispersion function 
into the negative imaginary frequency plane (r < 0) 



z-{f) 



7T 



1/2 



dx- 



x - f 



+ 2i^e~ f \ 9f(/)<0, 



and the plasma waves are covariantly described now for all complex frequencies by 



0"ll 



NL 



n q 3 fi 



Zl 



V2m 2 CU1U3 z 3 — Z\ 



Z1Z3 
Z3 ~ Zl 



J+(l 



where J is given by 



(46) 



(47) 



J = J + + 2%y/TT 



HI 



-^3) 



- r 2 

J3 



H( 



1-zl 



(48) 



A. Limit infinite speed of light c — > 00 

We finaly consider the plasma waves interaction in the formal limit of an infinitely large 
speed of light, c — > 00, which corresponds to the classical noncovariant nonrelativistic theory 
and where all waves have subluminal phase velocities. From Eq. (J44j) - (j48j) we obtain in this 
limit 



2^^(-^ 3 )e- (/ 3°°) 2 ) - f°° (Z + (A°°) + 2^^(-^i)e- (/r)2 )] +(1^2) (49) 



with 



v 2 _ 2^T 



m kiV th 
and which holds for all complex frequencies and for all wave-numbers. 



(50) 
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Assuming again |/°°| ^> 1, the exponential functions from the residues in (|49|) are very 
small (e - ^ 00 -* 2 <C 1) and the asymptotic approximation of plasma dispersion function in (}4"9~j) 
leads to 



nl,oo_. n oQ 3 N 3 (Nx + N 3 ) 



2m 2 c 



+ (!<-> 2), (51) 



which agrees with the standard results [|9j]] and where N = kc/ui is the index of refraction. 



VI. CONCLUSION 



On the basis of relativistic Vlasov-Maxwell equations we have developed a covariant 
kinetic formalism to determine the nonlinear plasma conductivity which allows us to find 
all the plasma waves nonlinearly excited in plasma. The general covariant form of nonlinear 
conductivity is provided first for any value of plasma temperature and for the whole complex 
frequency plane by a correct analytical continuation. Then, we have restricted the analysis 
to an appropriate relativistic particle distribution which is vanishing for particle speeds 
greater than speed of light. And in the limit of nonrelativistic plasma temperatures we 
have derived apparently for the first time the covariant nonlinear conductivity which is 
significantly different from the standard noncovariant nonrelativistic results. Only in the 
strictly unphysical formal limit of an infinitely large speed of light c — > 00 the covariant 
forms reduce to the standard expressions. 
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